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Abstract. We report on detailed Monte Carlo inversion analysis of the Lyman limit system observed at z a b s = 
2.917 in the VLT/UVES spectrum of the quasar HE 0940-1050. Metal absorption lines of carbon and silicon in 
three ionization stages and numerous atomic hydrogen lines have been analyzed simultaneously. It is found that 
in order to match the observations, the shape of the ultraviolet background ionizing spectrum of Haardt & Madau 
(1996) should be modified: a spectrum with a higher intensity of the emission feature at 3 Ryd is required. It is 
also found that synthetic galactic spectra (or different mixtures of them with power law spectra) cannot reproduce 
the observations, indicating that the stellar contribution to the ionizing background is negligible at z ~ 3. For the 
first time a very low carbon abundance of [C/H] = —2.93 ± 0.13 and the abundance ratio [Si/C] = 0.35 ± 0.15 
are directly measured in the Lyman limit system with iV(Hl) = 3.2 x 10 17 cm -2 . If the absorber at z a b s = 2.917 
provides an example of a pristine gas enriched by the nucleosynthetic products of early generations of stars, then 
the measured value of [Si/C] seems to indicate that the initial mass functions for these stellar populations are 
constrained to intermediate masses, M up < 25Mq. 

Key words. Cosmology: observations - Line: formation - Line: profiles - Galaxies: abundances - Quasars: absorp- 
tion lines - Quasars: individual: HE 0940-1050 



1. Introduction 

Metal abundance measurements in quasar intervening ab- 
sorbers provide insight into the chemical evolution of mat- 
ter over cosmological time scales. Since heavy elements 
were discovered in some Lya clouds at z ~ 3 by Cowie 
et al. (1995) and by Tytler et al. (1995), several mech- 
anisms explaining the origin of metals in the intergalac- 
tic medium (IGM) have been proposed. Nowadays it is 
supposed that at z > 10 the IGM was pre-enriched by 
the first Population III (Pop III) stars formed from gas 
with zero metallicity (e.g. Bromm et al. 2001; Nakamura 
& Umemura 2001) and enriched later due to disruption 
of low-mass protogalaxies (e.g. Madau et al. 2001) at 
6 < z < 10, and due to ejection of metals from massive 
galaxies at z < 6 (e.g. Aguirre et al. 2001; Scannapieco 
et al. 2002). It is clear that to see the signature of 
Pop III stars one must go to very low metallicities. The 
simulations of collapse and fragmentation of primordial 
gas clouds suggest that in the first generation the stars 
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may have been very massive, M* > 100M© (Abel et al. 
2000; Bromm et al. 2001). However, the calculations by 
Nakamura & Umemura (2001) produced a bimodal ini- 
tial mass function (IMF) for Pop III stars with a second 
peak at 1 — 2Mq. According to these calculations the first 
generation supernova (10M Q < M* < 35M©) and pair- 
instability supernova (140M Q < M* < 260M©) can pro- 
duce metallicity in the range 10~ 4 — 10 _3 Z©. Since the ele- 
ment yields and production rate differ significantly for the 
massive (10— 35M©) and very massive (140— 260A/©) stars 
(Woosley & Weaver 1995; Umeda & Nomoto 2001; Heger 
& Woosley 2002), the Pop III IMF can be constrained by 
measuring the abundance ratios in low metallicity cosmic 
objects. 

In the present study we report on the discovery of 
extremely metal-poor cloud at z abs = 2.917 toward the 
quasar HE 0940-1050 which reveals the lowest metallic- 
ity of [C/H]0 ~ —3 measured up to now at high red- 
shifts. All previously reported measurements of the metal 
abundances in the Lya clouds lie in the range [C/H] 
> -2.4 (Fan 1995; Songaila & Cowie 1996; Levshakov et 

1 Using the customary definition [X/H] = log (X/H) - 
log(X/H)0. Photospheric solar abundances are taken from 
Holweger (2001). 
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al. 2002a, hereafter LACM; Levshakov et al. 2003, here- 
after LAD WD). We note that the lowest carbon abun- 
dance measured in Galactic metal poor stars is [C/H] = 
-3.6 (Norris et al. 2001). 

While dealing with metallicity measurements in high 
rcdshift absorbers, one should take into account that the 
results are strongly dependent on several assumptions on 
physical parameters poorly known such as geometry of 
the cloud, gas density distribution, shape and intensity of 
the background ionizing radiation, ionizing mechanisms 
etc. Our results are obtained with the Monte Carlo inver- 
sion (MCI) algorithms described in detail in Levshakov 
et al.(2000, hereafter LAK), LACM, LADWD. Since this 
technique is relatively new, we briefly outline here its ba- 
sics. 

The main assumption is that all lines observed in an 
absorption system arise in a continuous absorbing gas slab 
of a thickness L with a fluctuating gas density and a ran- 
dom velocity field. Numerous cosmological hydrodynam- 
ical calculations performed in the previous decade have 
shown that the QSO absorption lines arise more likely 
in the smoothly fluctuating intcrgalactic medium in a 
network of sheets, filaments, and halos (e.g. Cen et al. 
1994; Miralda-Escude et al. 1996; Theuns et al. 1998). 
This is also supported by modern high resolution spectro- 
scopic observations: the increasing spectral resolution re- 
veals progressively more and more complex profiles. Thus, 
such an approach seems to be more physically justified as 
compared to the model of separate homogeneous clouds 
which the commonly used procedure is based on. This 
procedure consists of the deconvolution of the complex 
spectra into a set of separate Voigt profiles and the subse- 
quent estimation of ionization parameters (one for every 
subsystem) from the obtained ionic column densities. As 
shown in LAK and LACM, this simplified technique may 
produce incorrect metallicities exceeding in some cases 1 
dex. 

Further we assume that within the absorber the metal 
abundances are constant, the gas is in thermal and ioniz- 
ing equilibrium and it is optically thin for the ionizing UV 
radiation. The intensity and the shape of the background 
ionizing radiation is considered as an external parame- 
ter. Within the absorbing region the radial velocity v(x) 
and the total volumetric gas density nu{x) along the line 
of sight are considered as two continuous random func- 
tions which are represented by their sampled values at 
equally spaced intervals Ax. The computational proce- 
dure is based on the adaptive simulated annealing (see 
LAK, LACM), the fractional ionizations of different ele- 
ments are computed at every space coordinate x with the 
photoionization code CLOUDY (Ferland 1997). 

The MCI allows us to recover self-consistently the 
physical parameters of the intervening gas cloud; namely, 
the mean ionization parameter Uo, the total hydrogen col- 
umn density iVn, the line-of-sight velocity and density dis- 
persions (er v and a yi respectively) of the absorbing gas, 
and the chemical abundances Z a of all elements involved 
in the analysis. With these parameters we can further cal- 



culate the mean gas number density no, the column densi- 
ties for different species N a , the mean kinetic temperature 
Tkin, and the linear size L. Having these comprehensive 
information we are able to classify the absorber more reli- 
ably and hence to obtain important clues concerning the 
physical conditions in the intervening clouds. 

The structure of the paper is as follows. Sect. 2 de- 
scribes the data sets, the estimated parameters for the 
z a bs = 2.917 Lyman limit system are given in Sect. 3, the 
obtained results are discussed in Sect. 4, and our conclu- 
sions are reported in Sect. 5. 

2. Observations and data reduction 

Observations of the QSO HE 0940-1050 (B=16.6, z cm = 
3.06) from the Hamburg/ESO survey (Reimers et al. 1995; 
Rcimers et al. 1996; Reimers & Wisotzki 1997; Wisotzki et 
al. 2000) have been obtained with the Ultraviolet- Visual 
Echelle Spectrograph (UVES) on the Nasmyth focus of the 
ESO 8.2m KUEYEN telescope, second unit of the VLT at 
Paranal, Chile. 

The spectra were recorded with two different dichroic 
filters which allow to use the UVES blue and red arms si- 
multaneously as two independent spectrographs for both 
instrumental settings. One of the instrumental configu- 
rations covered the wavelength ranges 3300-3864 A and 
4790-6816 A with a wavelength gap between 5763 A and 
5846 A, while the other covered the spectral regions 3740- 
4983 A and 6715-10400 A with a wavelength gap between 
8530 A and 8677 A. 

Three exposures of 3600 seconds each one of 2700 sec- 
onds were obtained for each instrumental configuration 
over the nights 26-29 March and 3 April 2000, when the 
seeing, as given by the telescope guide probe, was between 
0.5 and 0.7 arcsec FWHM. 

The slit widths were set at 1 arcsec and the CCDs were 
read-out in 2x2 pixel binned mode, resulting in a spectral 
resolution between 6.6 and 7.1 km s _1 . This full width at 
half maximum of the instrumental profile was measured 
from the widths of emission lines of the Thorium- Argon 
lamp used for the wavelength calibration of the spectra. 

The data reduction was performed using the 
ECHELLE context routines implemented in the ESO 
MIDAS package. Flat-fielding, cosmic ray removal, sky 
subtraction, and wavelength calibration were performed 
on each spectrum separately. Typical r.m.s. of the wave- 
length calibration is < 1 mA. 

The observed wavelength scale of each spectrum was 
then transformed into vacuum, heliocentric wavelength 
scale. The single extracted spectra were then added to- 
gether using weights proportional to their S/N. Finally, 
the local continuum was determined in the average spec- 
trum by using a spline to connect smoothly the regions 
free from absorption features. The continuum for the Lya 
forest region was fitted by using the small regions deemed 
to be free of absorptions and by drawing an interpolating 
spline between them. 
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Fig. 1. UV-background ionizing continua used in the 
present calculations. The spectra have been normalized 
so that J v {Jiv = 1 Ryd) = 1.0. The solid curve shows the 
spectrum computed by Haardt & Madau (1996) at z = 3. 
The dashed curve is a galactic spectrum from Leitherer et 
al. (1999) shown in their Fig. 8(e) - spectral energy dis- 
tribution at 5 Myr with Z/Z Q = 0.001. The dotted line at 
3 Ryd shows the modification of the Haardt and Madau 
spectrum required by the present observations. The posi- 
tions of ionization thresholds of different Si and C ions are 
indicated by tick marks 

We confirm the identification of several metal 
absorption-line systems with z = 2.82, 2.32, 1.918, and 
1.06 performed from the low (4 A) resolution spectra of 
HE 0940-1050 by Reimers et al. (1995). These systems will 
be analyzed in detail in a separate paper. A new Lyman 
limit system (LLS) with z & bs — 2.917, found in the high 
resolution (~ 0.1 A) VLT/UVES spectra, is the subject 
of the present work. 

3. Lyman limit system at z a b s = 2.917 

The Lyman limit system at z a b s — 2.9171 toward the 
quasar HE 0940-1050 consists of hydrogen Lyman se- 
ries lines (from Lya up to Ly-20) and metal lines of C n 
A1334, Sin A1260 and A1526 (the Sin AA1190, 1193 lines 
are blended), Al n A1670, Cm A977, Sim A1206, Nm 
A989, Civ A1548 and A1550, Siiv A1393 and A1402. The 
doublet Ovi AA1031,1037 is hopelessly blended and no 
conclusions about this important element can be made. 
Additionally the quasar spectrum beyond the Lyman limit 
(A bs < 3572.4 A) shows residual intensity I\ ~ 0.14 giv- 
ing us a model independent estimation of the total neutral 
hydrogen column density of about 3 x 10 17 cm~ 2 . 

The first attempts to inverse the system under study 
were made with the HM ionizing spectrum correspond- 
ing to z = 3. In spite of many trials with different initial 
conditions it turned out to be impossible to gain a good 
fitting for all available lines: the Sin A1260 line was sys- 



Table 1. Physical parameters of the z a b s = 2.917 Lyman 
limit system toward HE 0940-1050 derived by the MCI 
procedure 
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a Z x = X/H, [Z x ] = log(Zx) - log(Zx)©. 

6 A1 solar abundance from Grevesse & Sauval (1998). 

Internal error is 15%. 

internal error is 20% 



tematically overestimated, whereas the Siiv A1393, 1402 
doublet - underestimated. This was the reason why we 
started to modify the shape of the background ioniz- 
ing spectrum. After numerous experiments with different 
power law spectra, galactic spectra from Leitherer et al. 
(1999) and combinations of them it was found empirically 
that the self-consistent fitting of all observed lines can be 
achieved if we adopt the HM spectrum with about 4 times 
enhanced intensity between 2 and 3 Ryd (shown by the 
dotted line in Fig. 1). Can this emission peak be present in 
the real intergalactic UV background ? Yes, indeed. The 
bump at 3 Ryd in the HM spectrum is due to He II Lya re- 
combination radiation produced by the photoionized Lya 
forest clouds. The primary quasar UV continuum spectra 
were adopted by HM in the form of power law v~ 15 (for 
A < 1216 A). But it is known that the far UV region of 
quasar radiation sometimes contains rather strong emis- 
sion lines. For instance, quasar He II Lya emission at 3 Ryd 
was observed in spectra of HS 1700+6416 (Davidsen et al. 
1996; Reimers et al. 1998), Q 0302-003 (Jakobsen et al. 
1994; Heap et al. 2000), and in the composite QSO spec- 
trum (Telfer et al. 2002), the last spectrum also shows ad- 
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Fig. 2. Hydrogen and metal absorption lines associated with the z a b s = 2.917 LLS toward HE 0940-1050 (normalized 
intensities are shown by dots with la error bars). The zero radial velocity is fixed at z = 2.9171. Smooth lines are 
the synthetic spectra convolved with the corresponding point-spread spectrograph function and computed with the 
physical parameters listed in Table 1. Bold horizontal lines mark pixels included in the optimization procedure. The 
normalized Xmin = 0-97 (the number of degrees of freedom v = 553). Profiles of Ly-/3 (Av < — 70 km s -1 ), Ly-7 
(Av < - 80 km s- 1 and Av > 40 km s" 1 ), Ly-7 (Av > 35 km s" 1 ), Ly-15 [Av < - 50 km s" 1 ), and Cm A977 
(Av < — 35 km s _1 ) are contaminated by other absorption lines (presumably from the Ly-a forest). The blue wing of 
Ly-6 (Av < — 60 km s _1 ) is blended by the Ly-/? line from the z abs = 2.5536 system. The red wing of Ly-15 (Av > 20 
km s _1 ) is blended by the Ly-/3 and Ly-5 lines from, respectively, the 2 a b s = 2.4962 and 2.8237 systems (the systems 
with z a b s = 2.4962 and 2.5536 are identified in I. I. Agafonova et al. 2003, in preparation) 



ditional three wide emission features between 304 A and 
500 A. Thus the strong bump at 3 Ryd can be explained 
qualitatively. However, to estimate its exact shape and 
the amplitude, calculations accounting for the presence of 
emission lines in far UV QSO spectra and for radiative 



transfer effects in an inhomogeneous universe are needed. 
From the analysis of the z a b s = 2.917 system we only con- 
clude that to match all Si and C lines, about 4 times more 
photons with energies between 2 and 3 Ryd are required 
(the dotted line shown in Fig. 1 is only one possible ex- 
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Fig. 3. Portion of the VLT/UVES spectrum of HE 0940-1050 in the region of the Lyman limit absorption at z a bs — 
2.917. The panel illustrates the higher order hydrogen Lyman series lines from the z a ba — 2.917 LLS. Smooth curve 
represents a synthetic spectrum calculated with the physical parameters listed in Table 1. Continuum depletion between 
3538 A and 3555 A is due to the Lya absorption from the damped Lya system at z a b s = 1.918 



ample, the emission bump at 3 Ryd may be broader and 
lower as well). 

Due to the large number of ionic transitions in the 
z a b s = 2.917 LLS and due to the very high S/N ratio this 
system is a sensitive probe of the shape of the background 
photoionizing radiation. For instance, it was impossible to 
describe self-consistently all observed profiles of hydrogen 
and metal absorption lines using any of the composite UV 
spectra (A < 912 A) of young galaxies from Leitherer et al. 
(1999) (a typical galactic spectrum is shown in Fig. 1 by 
the dashed line) or any of the combined spectra calculated 
by mixing (with different weights) a spectrum of the young 
galaxy with the HM UV background. From this, we do 
not confirm the result by Steidel et al. (2001) that the 
Lyman-break galaxies produce a dominated photoionizing 
background at z <~ 3. 

The results obtained with the MCI and the modified 
HM spectrum (Fig. 1) are presented in Table 1 and illus- 
trated in Figs. 2-4 (for details of calculations see LAK and 
LACM) . Parts of line profiles included in the least-squares 
minimization are marked by horizontal lines in each panel 
in Fig. 2. We also tried to probe the presence of deuterium 
absorption but found that the wide Lya profile is not sen- 
sitive to any additional D absorption. The blue wing of 
the Ly/3 line is strongly blended and this line cannot be 
used in the D/H estimations. 

Fig. 2 shows that the observed profiles (portions free 
from blending) are well represented. A self-consistent fit- 
ting of all available lines has revealed that the blue wing 
of the C in A977 line is partly contaminated by some Lya 
forest absorption. 

The profiles of the Alii A 1670 and Nm A989 were cal- 
culated in a second round using the velocity v(x) and den- 
sity nn(x) distributions already obtained and the metal- 
licities chosen in such a way that the synthetic spectra did 
not exceed 1 <r deviations from the observed intensities. 



The Lyman series lines higher than Ly-12 are lying on 
the right wing of the damped Lya (DLA) identified at 
z a bs — 1-918 (Fig. 3). Using the neutral hydrogen column 
density in the z a b s = 2.917 LLS, we estimated the neutral 
hydrogen column density for the DLA system of iV(Hl) 
~ 1.0 x 10 20 cm~ 2 (this DLA system is described in detail 
in M. Centurion et al. 2003, in preparation). 

Fig. 4 illustrates the distributions of the gas density 
and the radial velocity inside the z a b s = 2.917 LLS. It is 
clearly seen that the co-existence of lowly and highly ion- 
ized species in this system is due to density fluctuations. 
A wide range between x = 0.15 and 0.37 with very low 
gas density is also responsible for the blue-side asymmetry 
seen in the metal line profiles in Fig. 2. 

4. Discussion 

According to the recovered parameters, the LLS at z a bs = 
2.917 is a very large could (L > 100 kpc) with very low 
metal content. What could it be related to ? 

In our code L is calculated from the ratio L = Na/no, 
and no oc (1 + <r 2 ) J912/U0 (here J912 is the Lyman-limit 
specific flux in erg cm~ 2 s _1 Hz -1 sr _1 ). Thus the inten- 
sity of the background radiation can influence the mean 
gas density and, hence, the absorber size. We used in our 
calculations J912 = 0.4 x 10~ 21 erg cm~ 2 s _1 Hz -1 sr _1 
given by Haardt & Madau (1996) for the mean background 
flux at z = 3.0. This value is poorly known and may be 
affected by local sources. But since we have no clues on 
how this value might be changed for this particular case, 
we will rely on the estimated linear size. 

For linear sizes as large as hundreds of kpc the velocity 
differences due to the Hubble expansion, A«h, become 
already noticeable. At z a b s = 2.917 and L = 140 kpc the 
velocity difference is Avh = 81 km s _1 , assuming Q = 1 
and the present day Hubble constant H = 75 km s _1 
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Fig. 4. Computed velocity (upper panel) and gas density 
(lower panel) distributions along the line of sight within 
the z abs = 2.917 absorber toward HE 0940-1050. Shown 
are patterns rearranged according to the principle of min- 
imum entropy production rate (see LACM) 



tively, in the Pair Instability SNe modelQ. Similar results 
are reported in Heger & Woosley (2002) who also pre- 
dicted a pronounced deficit of Al Compared to Si ([Al/Si] 
< — 1.5) in the mass range 140 — 260M Q . 

In our system we measured [Si/C] = 0.35 ± 0.15 and 
[Al/Si] = —0.18 ±0.13. These values certainly rule out the 
enrichment by very massive stars explosions and restrict 
the stellar masses by ~ 25M Q . It should be noted that the 
abundance ratio [Si/C] ~ 0.3 — 0.4 is in line with other 
measurements in high redshift absorbers (e.g., Songaila 
& Cowie 1996; Levshakov et al. 2002b, 2003). If, never- 
theless, the Pop III stars were massive (M > 140M Q ), 
then the mctallicity produced by their explosions should 
be much lower than 10~ 3 Zq. 

We also estimated an upper limit on [N/C] < —0.4. 
Nitrogen is usually assumed to be mostly produced by 
the intermediate mass stars (4 < M/Mq < 8) whereas 
carbon production is dominated by the massive (M/Mq > 
8) stars (Henry et al. 2000). Unfortunately, the obtained 
upper limit is too high to allow us to conclude about the 
possible nature of nitrogen in the z a hs = 2.917 LLS: the 
Niii A989 line lies in the Lya forest where the signal-to- 
noise ratio is only about 25. 



Mpc^ 1 . The corresponding velocity dispersion caused by 
the Hubble flow is then oh = Ai>h/\/I2 = 23 km s _1 . 
The velocity dispersion of the bulk material deduced by 
the MCI is er v = 50 km s _1 , implying that the line-of-sight 
peculiar velocity dispersion is about 44 km s _1 . 

Large size, small peculiar velocity dispersion, and very 
low metallicity may indicate that the z a b s = 2.917 absorber 
is hosted by either an external galactic halo or by some 
large-scale structure object like a filament, and may con- 
sist of the gas pre-enriched by the first generations of stars. 
The sizes of the galactic halos or the filaments have been 
recently directly measured at z < 1. For instance, ob- 
servations of Civ absorption lines by Chen et al. (2001) 
revealed that all these lines can be associated with galax- 
ies seen up to 180 kpc away. Penton et al. (2002) reported 
on large size Lya absorbers aligned along a possible fila- 
mentary structure of ~ 20h^Q Mpc long and ~ lftfo 1 Mpc 
wide. 

As already mentioned above, the metal content mea- 
sured in this LLS may be considered as primordial, i.e. 
produced by the explosions of early generations of stars. 
Then we may use the relative metal abundances to restrict 
masses of these hypothetical Pop III objects. Theory pre- 
dicts that the abundance ratio [Si/C] rises to highly su- 
persolar values with increasing masses of the Pop III stars. 
The calculations of the nucleosynthesis patterns of metal- 
free stars with various masses (e.g. Umeda & Nomoto 
2001) show that [Si/C] = 0.34, 0.42, 0.70 for, respec- 
tively, M/Mq = 13, 15,20 (Type II SNe); [Si/C] = 0.41 - 
0.76 for the Pop III 25M© model with explosion energies 
E exp = 10 51 — 10 52 erg (Fe Core Collapse Hypernovae); and 
[Si/C] = 1.49 and 1.70 for M/Mq = 170 and 200, respec- 



5. Summary 

We have analyzed the VLT/UVES spectrum of the quasar 
HE 0940-1050 and deduced the physical properties of the 
Lyman limit system at z a hs = 2.917. The main conclusions 
are as follows : 

1. For the first time, a very low carbon abundance of 
[C/H] = -2.93 ± 0.13 is directly measured in the LLS 
providing an example of a gas cloud probably enriched 
by the early generations (z > 3) of stars. 

2. The analyzed Lyman limit system has the line-of-sight 
size L of about 140 kpc. Its velocity dispersion is only 
about 50 km s _1 . This system may be hosted by an 
intergalactic filament structure or an external galactic 
halo. 

3. The measured abundance ratios [Si/C] = 0.35 ± 0.15 
and [Al/Si] = —0.18 ± 0.13 are in line with the as- 
sumption that the initial mass functions for early gen- 
erations of stars are restricted to intermediate stellar 
masses, probably M up < 25M Q . The relative abun- 
dances of [Si/C] and [Al/Si] do not reveal any signa- 
ture of very massive stars from the early generation. 

4. The estimated shape of the photoionizing background 
radiation does not confirm that the Lyman-continuum 
radiation escaping from the young galaxies produces 5 
times more H-ionizing photons per unit comoving vol- 
ume than QSOs at z ~ 3 as suggested by Steidel et 
al. (2001). None of the far UV galactic spectra calcu- 
lated by Leitherer et al. (1999) can produce profiles of 

2 The [Si/C] ratios from Umeda & Nomoto (2001) which are 
normalized to the solar abundance system given by Anders & 
Grevesse (1989) are adjusted to the new solar photospheric 
values from Holweger (2001). 
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carbon and silicon ions observed in the LLS at z a b s = 
2.917. 
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